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A family of the pair hopping models exhibiting the incompressible quantum liquid at fractional 
filling l/m D is constructed in D dimensional lattice. Except in one dimension, the lattice is the 
generalized edge-shared triangular lattice, for example the triangular lattice in two dimensions and 
tetrahedral lattice in three dimensions. They obey the new symmetry, conservation of the center- 
of-mass position proposed by Seidel et al.. 1 The uniqueness of the ground state is proved rigorously 
in the open boundary condition. The finiteness of the excitation energy is calculated by the single 
mode approximation. 



I. INTRODUCTION 

The featureless Mott insulator blocks the charge trans- 
portation due to the strong electron-electron interaction 
and the ground state exhibits no symmetry breaking. 
Their existence is very rare both in the experimental sys- 
tems and in the theoretical models. There are two exam- 
ples in the two dimensional systems. One is the famous 
fractional quantum Hall effect, where the ground state 
is an incompressible liquid when the filling factor is 1/q 
(q = odd integer). The other example is the quantum 
dimer model in the two dimensional triangular lattice, 
where the ground state is a disordered dimer liquid that 
separates the excitations by a finite energy gap. 

Very interestingly, Seidel et al. pointed out that 
these two systems actually belong to the same type 
of Hamiltonian. 1,2 Namely, their Hamiltonians preserve 
both the center-of-mass momentum and the center-of- 
mass position. They showed that the fractional quantum 
Hall system in the lowest Landau level on the torus de- 
scribed by the following pseudo-potential Hamiltonian 3 

H = J d 2 rd 2 r'V 2 <5(r-r')^ + (r)^ + (r')V'(r')V'(i-) (1) 

can map to a pair hopping model in one dimensional 
lattice 

H= ]T f*(x)f(y)C R+x C R _ x C R - V C R+V 

R,x,y 

f{x) =(t 3/2 ^(i-nI)e- K2(l - ni » 2 , (2) 

n 

where 

L = L x L y /2Trl%, k = 2-Kle/Ly, (3) 

L x and L y are the linear dimensions of the torus, and 
Ib = \JhcjeB is the magnetic length. It can be obvi- 
ously seen that Eq.(2) describes the hopping that pre- 
serves the center-of-mass position. Two electrons anni- 
hilated at R ± y hop to R ± x with the effective hopping 
range 1/k and center-of-mass position preserves at R. 
Due to this beautiful property, the ground states of the 
1/q fractional quantum Hall liquid with the g-fold de- 
generacy can be labeled by the q different center-of-mass 



positions. Furthermore, Seidel et al. also showed that 
the ground state of Eq.(2) is a charge density wave with 
the amplitude ~ e~ c / K where c is a constant of 0(1), and 
the energy gap is finite for any finite k. For very small 
k, where Eq.(2) becomes the long range hopping model, 
the charge density wave amplitude is exponentially small. 
One can safely say that the ground state describes a fea- 
tureless Mott insulator without local order parameter. 

The new symmetry of the center-of-mass position con- 
servation paves a new way to search the featureless Mott 
insulator. Recently, inspired by the higher dimensional 
generalization of the quantum Hall effect 4 and the Hal- 
dane's pseudopotcntial Hamiltonian, Chern et al. con- 
structed a model for the incompressible liquid in the two- 
dimensional triangular lattice. 6 They showed that the 
ground state is unique without local order parameter. 
They also computed the excitation gap using the sin- 
gle mode approximation. This paper will serve as the 
extended version of that Letter. For this purpose, we or- 
ganize the paper in the following: In the section II, we 
review the Haldanc construction of the pseudopotcntial 
method on the two sphere. We will show explicitly that 
the Haldane pseudopotential Hamiltonian can map to a 
long-range pair hopping model with the conservation of 
the center of mass position. In section III, we provide 
the detail calculation of our previous Letter and gener- 
alize it to the SU(N) model. In section IV, we provide 
the detail calculation of the single mode approximation 
in the SU(3) case. Finally, we conclude and summarize 
in the section V. We also include several appendixes for 
the readers to follow the group-theoretical method easily. 

II. ONE-DIMENSIONAL LATTICE MODEL FOR 
THE INCOMPRESSIBLE QUANTUM LIQUID 

The quantum Hall effect can be considered on the 
two sphere subject to the uniform magnetic field by the 
U(l) magnetic monopole at the center. 5 In this case, the 
Laughlin wave function becomes the exact ground state 
of the Eq.(l). We will show that it can map to a pair hop- 
ping model with the conservation of the center-of-mass 
position and long range hopping integral. 
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In the presence of the U(l) magnetic monopole flux, 
the single-particle wavcfunction is described by the 
monopole vector spherical harmonics 7,8 which can be de- 
noted by the SU(2) |Z, m > state, where I can be integers 
or half-integers and m is the magnetic quantum number. 
The Landau level spectrum is given by 5 



Eu 



2MB? 



1(1 + 1) 



(4) 



where M is the mass of the electrons and R is the radius 
of the sphere. If the total magnetic flux is 25, I = S + k, 
where k is the Landau level index. Because k is only an 
integer, 5 is either integer or half-integer. Each Landau 
level has 2(5 + k) + 1 degeneracy. In the lowest Lan- 
dau level, k = 0, the single-particle wavefunction can be 
written as 



(25)! 



{S + i)\(S-i)\ 



u s+ *v s - 



(5) 



where i = —5, —5 + 1, .., 5 — 1,5 and the (u, v) is the 
two-dimensional complex spinor given by 



| e -i*/2 



(6) 



where 8 and <p parameterizing the sphere are known as 
the polar and azimuthal coordinates. This is the special 
property of the 111 that the single-particle wavefunction 
can be completely described only by one quantum num- 
ber i. Therefore, one can treat the configuration space 
as the ID chain with the number of site 25 + 1 and the 
lattice site is labelled by i. If the filling factor v — 1, the 
number of the particles N = 25+1. Then, the many- 
body wavefunction is the Slater determinant given by 



JV 



(7) 



k<l 



where k and I are the particle indices. It is easy to see 
that Eq.(7) is the unique many-body fermionic wavefunc- 
tion for v = 1. 

The quantum Hall state with filling factor v = 1/m 
celebrated as the Laughlin wavefunction 9 can be written 



as" 



(8) 



k<i 



where m is an odd integer. In Eq.(8), the maximum 
power of u k becomes 2m5, indicating that the single- 
particle wavefunction given by the Eq.(5) is in the spin- 
mS state and therefore the dimension of the single- par- 
ticle states (or say the lattice) is 2m5+ 1. While keeping 
the number of particles = 25 + 1 the same, Eq.(8) 
describes a state with filling factor 



25+1 
2m5 + 1 ' 



(9) 



which goes to 1/m as 5 goes to infinity. Because 5 scales 
as R 2 shown in the Eq.(4), infinite 5 indicates to take the 
thermodynamic limit. 

Eq.(8) is the unique ground state of the following 
many-body Hamiltonian 



Hi 



1 



k<m-2 
(ij) k—1, odd 



E «*r< 



2mS-k 
ij 



(10) 



where a k are positive-definite and P+™ are the pro- 
jection operators that project two spin-m5 states to the 
two-body states of total spin 2m5 — k for the pair (ij). 
We note that k can be only the odd integers because the 
two-body 2m5 — k states with odd k are antisymmetric 
upon particle exchanges. In Eq.(8), the term with the 
maximum power of u^Uj for any pair (ij) is 



( Ui Uj) 



m(2S-l) 



(UiVj 



(11) 



which indicates that no two-body 2m5 — k states for k < 
m — 2 for any pair (ij). Therefore, Eq.(8) is the zero 
energy state of Eq.(10). On the other hand, because 
there is no two-body 2m5 — k state for k < m — 2 for any 
pair (ij), any ground state wavefunction \ must have the 
following form 

N 

X ~ -;u N ,v N ) Y[(u k vi - uiv k ) m , (12) 

k<l 

where / is the symmetric function for any pair exchange. 
For the single-particle wavefunction to be described by 
I = mS, the power of u k for each particle in Eq.(12) has 
to be 2m5. However, the factor 



N 



Y[(ukVi - uiv k ) r 



(13) 



k<l 



in Eq.(12) already exhausts the quota of the power of u k - 
Then, / can only be a constant. \ ~ *m an d Eq.(8) is 
indeed the unique ground state of Eq.(10). 

The SU(2) spin model of Eq.(10) can be formulated as 
the lattice hopping model. As we mentioned, the single- 
particle state can be labelled by one quantum number i. 
Let us denote 



\mS,i) ->ct|0) 
Then, Eq.(10) can be written as 



(14) 



Ul = a EE fc ' 9) C l C I-fc C «-9 C 9- (I 5 ) 

n k,q 

for m — 3, where F(n, k, q) can be obtained by the SU(2) 
Clebsh-Gordan coefficient 

F(n,k,q) = f*(n,q)f(n,k), 

M (-l) 2n y35(2fc-n)Vl25-l 
f(n,k) = 



^(35 - k)(3S + k)(3S + k- n) 
7(65-71-1)1(65 + 71-1)! 
7(35 -k + n)(125 - 1)!(35 - k - 1)! 

7(65-1)! 

7(35+fc-l)!(35+A:-n-l)!(35-A;+n-l)! 



(16) 
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It can be easily seen that Eq.(15) is a center-of-mass 
conserving hopping. Two electrons annihilated with the 
center-of-mass position |(fc + n — k) = | will be created 
in pair with the same center-of mass position. The rea- 
son that the hopping matrix element F is n-depcndent is 
due to the presence of the open boundary. Although F 
has long range hopping, the hopping matrix elements is 
exponentially small as the relative distance is comparable 
to the system size. 

To analysis hopping range, let us consider the hopping 
with the center-of-mass to be at the origin, namely n = 0. 
Then /(n, k) = g(k), where 



,(*) = ^kVssiies-iy.r 

(fc 2 -9S 2 V(12S-2)!(3S-fc-l)!(3S+fc-l}! 

g(k = 0) = indicates that no two electrons can be 
created (annihilated) at the same sites. For S = 30, g(k) 
is calculated in the Fig.(l). Taking k = 3S, 



g(3S) ~ VSe- 6810 * 2 , 



(18) 



which has an exponentially small tail. A proper definition 
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FIG. 1. The g{k) versus k for S = 30. 

of the hopping range can be defined as the half width of 
the hump of the \g(k)\ 2 . The result that the hopping 
range scales with S is given in the Fig. (2), where we 
show the log-log relation between the hopping range and 
the size of the system L = 6S + 1. The trend line in 
the Fig. (2) is the best fit given by |logL — 0.9, which 

indicates that the hopping range scales as L 1 / 2 . 

In summary, using the SU(2) weight space, we have 
considered the 2D fractional quantum Hall effect as a 
ID lattice model with the open boundary condition. 
The weight number, labeling the lattice site, is the 
^-component coordinate of the coherent state on the 
sphere. That the fractional quantum Hall state does not 
exhibit any density- wave order on the sphere implies that 
the ground state has no long range order in the ID lattice 
model. On the other hand, the sphere has no edge. The 




FIG. 2. A Log-Log plot of the width in \g(k)\ 2 as 
a function of L = 6S + 1. We have calculated for 
S=0,l,2,.., 10,20,30,. .,100,200,300,.. ,600. 



energy spectrum of the Eq.(15) does not include the edge 
modes that occur in the open 2D plane system. Thus, 
there is no gapless excitation corresponding in this model. 
Note that although the fractional quantum Hall state is 
the incompressible liquid shown in the experiments, theo- 
retically the finiteness of the excitation gap is only proved 
by the single mode approximation 10,11 . 



III. HIGHER DIMENSIONAL LATTICE MODEL 
FOR THE INCOMPRESSIBLE QUANTUM 
LIQUID 

There is a natural generalization of Eq.(5) to the 
SU(N) coherent state, which is labelled by (N — 1) quan- 
tum numbers without degeneracy. The SU(N) coherent 
state is known as the SU(N) (p, 0, .., 0) multiplet. Similar 
to the SU(2) case, we will construct a family of the faith- 
ful lattice models for the incompressible quantum liquid 
using the SU(N) coherent states. We shall start with the 
SU(3) for the pedagogical purpose. 

The SU(3) coherent states are given by 



mi!m2!iTi3! 



u mi v m2 w m3 , m 1 +m 2 +m 3 =p (19) 



which forms the multiplet described the SU(3) (p, 0) rep- 
resentation and m,i are integers. The (u, v, w) in Eq.(19) 
is the complex spinor which can be represented by 



tp a (zi,z 2 ) 




1 



y/l + Z\Z\ + Z 2 Z 2 



(20) 



where Zi are the complex numbers to parameterize CP 2 
and a is the spinor index. As shown in the Appendix, 
the SU(3) (p, 0) states is the lowest Landau level in the 
quantum Hall problem in CP 2 . The SU(3) Cartan sub- 
algebra contains two generators T 3 and T 8 . The states of 
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Eq.(19) acquire the coordinates in the T3 and Tg, space, 
which form the two-dimensional triangular lattice shown 
in the Fig. 3. The number of sites of the two-dimensional 



OOO QOOO 



OOO 



O O 



OOO 



O O 



OOO 



O 



9 



where the operator p^ mp ~ 2q ' q } operates on the direct 
product states of two spins i and j and projects them onto 
the (2mp — 2q, q) states, and n q are positive-definite. The 
SU(3) two-spin states are the generalization of the SU(2) 
angular momentum addition. The direct product of two 
SU(3) multiplets can be also block-diagonalized such that 
each block is described a SU(3) multiplet denoted by two 
integers (a, b). In this case, the direct product of two 
SU(3) (mp, 0) multiplets can computed as 

mp 

{mp,0) x (mp,0)| a = (2mp-2q,q) (26) 

q— l,odd 

where the subscript "a" denotes the antisymmetric com- 
bination and k are the odd integers. Furthermore, ev- 
ery SU(3) multiplet can be also block-diagonalized by its 
SU(2) subgroup. In our case, the SU(3) (p, 0) multiplet 
can be decomposed as 



FIG. 3. The weight space of (6,0). T 3 and T 8 are the Gel- 
Mann matrices forming the Cartan subalgebra. 

lattice representing the SU(3) (p, 0) multiplet is given by 
1 



d(p) = -(p+l)( P + 2) 



(21) 



If the number of electrons N = d(p), the many-body 
fcrmionic wavefunction is the Slater determinant given 

by 



* = 



p-1. 

2 



u 



V2 



p v— 1 

u% u y N V N 



id, 



w 



N 



(22) 



up to the normalization constant. It is also not hard 
to see that Eq.(22) is the unique fcrmionic many-body 
wavefunction when TV = d(p). Next, let us consider the 
natural generalization of Eq.(8) written by 



'■N 



L 

V2 



p-1 



Wo 



'N 



(23) 



Because the highest power of u for each particle in 
Eq.(23) is mp, the new coherent state is described by 
the SU(3) (mp,0) multiplet. Therefore, Eq.(23) is the 
many-body state with the filling factor 



(p+l)(p + 2) 
(mp + l)(mp + 2) 



(24) 



which becomes 1/m 2 in the thermodynamic limit. We 
propose the following SU(3) spin Hamiltonian 



q<m—2 
9=1 > odd 



K P 

K q ^ij 



(2mp—2q.q) 



(25) 



(P,O) = 



(27) 



fc=0 



which can be easily checked by the counting the dimen- 
sionality. 

The Eq.(23) is the zero-energy state of Eq.(25). We 
first look at the two-spin state of Eq.(22). Because 
all sites are occupied, the SU(3) two-spin state with 
the maximum SU(2) quantum number is the (2p — 2, 1) 
multiplet. Because Eq.(23) is the (to- fold) product of 
the Eq.(22), the SU(3) multiplet with the maximum 
SU(2) quantum number that the direct product of to 
(2p — 2, 1) multiplets can yield is (2mp — 2m, to). There- 
fore, due to the SU(3) symmetry, the two-spin states in 
Eq.(23) does not contain the multiplets (2mp— 2k, k) for 
k < to — 2. The Hamiltonian is the SU(3) generaliza- 
tion of the Eq.(10). Our generalized Laughlin wavefunc- 
tion of Eq.(23) is the zero-energy state of the generalized 
Hamiltonian of Eq.(25). Moreover, it is also the unique 
ground state. In the Letter 6 , we have demonstrated rig- 
orously that Eq.(23) is the unique ground state of the 
Hamiltonian of Eq.(25). Some supplemental details will 
be discussed in the Appendix. 

Eq.(25) can be written as a lattice hopping model. 
Similar to Eq.(14), the electron creation operator can be 
defined as 



\(mp,0);j,j 3 



c L 3 )i°>' 



(28) 



where the states are denoted by the quantum numbers of 
the SU(2) subgroup. This coordinate system is equiva- 
lent to the quantum numbers of T 3 and Tg, since there is 
no degeneracy on the lattice. Using Eq.(28), the Hamil- 
tonian for to = 3 can be expressed by 



j,L,L$ 1,1$ k,k3 



L,L 3 
3 



c (i,h) c U-i+h,L 3 -i 3 fU-k+h^-k 3 ) c (k,k 3 )' ( 29 ) 

where n is a positive number and F?', L ' L3 is the SU(3) 
Clebsh-Gordan coefficient from two (3p, 0) multiplets to 
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the (6p — 2,1) subspace shown in Eq.(26). In Eq.(29) 
the center-of-mass position is conserved at (j ' + 1/2, L3] 
in the pair-annihilation and pair-creation process. 

The Hamiltonian of Eq.(29) also describes a long- 
ranged hopping process. The i^'j ' 3 can be computec 
exactly as 12,13 



F; 



jt L,L 3 _ (3p)! /(6p-l)2!(6p+i-i-2/-i) 



i.h 



(6p)! 



(3p-2j + 2Z-l) 



(Qp-L-j- §)!(2j + l)! V(2j - 2? + 2)!(2Z + 1)! 
(3p-2l)! 



(-l)3+J+i3-2Z^2L + 1 
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where t = j — I and the Wigner 6-j and 9-j symbols are 
explicitly used. In Fig. 4, we plot the F/f' 1 " 3 for p = 300, 



j = 3p - |, L = 3p - 1, and Z = f 
of Z3. It is the pair-hopping integral with the relative 
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FIG. 5. A log-log plot of the hopping range to p from 30 to 
600. The straight line is the best fit. The vertical axis is the 
log of the hopping range and the horizontal one is log p. The 
hopping range scales as p5 
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p,3p- 



,3p - 1 '°,p = 300, i 3 = l..* 



4Z§ . Because two electrons are not on the 



distance y^f 

same rows, Ffc : 3 is not zero at I3 — 0. The hopping 
range can be also defined as the half- width of IF 1 ! 2 . In 
Fig. 5, we show the log-log relation between the hopping 
range and p. The result suggests that the hopping range 
scales as ^fp. Similar to the SU(2) case, when the relative 
distance is the same order of the system size, the pair- 
hopping integral is exponentially decayed as 



F 



,p,3p-f ,3p-l,0 

3p 3p 



\/6p(6p - 1) e 



3p log 2 



(31) 



In the later section, we will show the existence of the 
finite excitation gap within the single-mode approxima- 
tion. In the presence of the energy gap, the uniqueness 
of the ground state implies that it is an incompressible 
quantum liquid. Here we shall consider the more general 



Our current formalism can be generalized to the SU(N) 
case very easily. The SU(N) fundamental spinor given by 



f Ul \ 

u 2 



V u N ) 
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Z-2 



\ZN J 



(32) 



can be used to parameterized CP^ 1 , where Zi are com- 
plex numbers. The SU(N) coherent states given by 



PI 



nf =1 



N 
t=l 



m l =p 



(33) 



are described by the SU(N) (p,0,0, ..,0) multiplet. A 
general SU(N) multiplet is labelled by N — 1 integers 
(m,..,n N ). The rank of the SU(N) group is N-l. There- 
fore, Eq.(33) represents a (A^— l)-dimensional lattice. As 
we have seen, it is a triangular lattice in two dimensions, 
and in three dimensions it is a tetrahedral lattice, and so 
on. One of the common features is that they are all frus- 
trated lattices. The number of the lattice is also related 
to p, namely the maximum power of the Uj, given by the 
following formula 



d N (p) = 



1 



(N-l)\ 



(p+\){p + 2)-ip + N-\). (34) 



When the number of the electrons is equal to the number 
of sites, A^ = dpf(p), the fermionic many-body wavefunc- 



6 



tion is given by 



N 
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(35) 



where Ujj is the i component of spinor in Eq.(32) for the 
j th electron. Similarly, let us consider the wavefunction 
^J}. Because the maximum power of u\ for each electron 
becomes rap, it describes a state with fractional filling 
factor 



v N = 



0+l)(p + 2) ■ jp + N- 1) 
(mp + l)(mp + 2) • -(mp + N - 1) 



(36) 



which approaches to I/to^ -1 or l/m D , where D is the 
dimensionality, in the thermodynamic limit. 

The Hamiltonian for the is the zero-energy state 
is given by 



q<m-2 



H=\Y. E * 



p. 



(2mp-2q,qfi,0,..fi) N 



(37) 



i^j q— 1, odd 



where n q are positive and p^ mp ~ 2q ' q ' ' '--'°^ N projects two 
SU(N) (mp,0, ..,0)jv spins to the (2mp-2q, q, 0, 0, ..,0)jv 
spin. The argument for to be the zero-energy state 
is quite similar to the SU(3) case. We also leave the 
proof for the non-degeneracy of the ground state in the 
Appendix. 

Each the SU(N) coherent state is non-degenerate. 
States in the SU(N) (p, 0, .., 0)jv multiplet can be labelled 
by its SU(N-l) subgroup following the relation 



(p,O,O,..,OV = 0(a:,O,..,O)j V _i 



(38) 



x=0 



From Eq.(27) and Eq.(38), there is an inductive relation 
to label the states in the SU(N) (p, 0, ..,0)jv multiplet. 
In other words, they can be labelled by the SU(N-l), 
SU(N-2), .., and SU(2) subgroups namely 



|(mp,0,0, ..,0) A r;a; 1 ,X2,..,a; A r_ 1 > 



(39) 



where Xi denote the SU(i) (xi,0, ..,0)j multiplets. Using 
this coordinate system, the electron creation operator is 
defined by 



|(p,0,0,..,0)iv;a;i,a;2,--,a;jv- 



Let us denote x — {x\,X2, ..,Xn-i). Then, similar to the 
Eq.(25), Eq.(37) can be written as 



H= 



K 



E 

ReQmp-2,1,0, ..,0) N V> X 



J2F(R,y)F(R,x)c^,c^c s (41) 



for m = 3, where F(R 1 x) is the Clebsh-Gordan co- 
efficient from two SU(N) (3p, 0, 0, .., 0)jv spins to the 
(6p — 2, 1,0, ..,0)jv spin. The only tricky point in 



Eq.(41) is the center-of-mass R. Because the SU(N) 
(6p — 2, 1, 0, 0, .., 0)jv multiplet contains degeneracy in its 
states, the (N — l)-dimensional vector R is the function 
of the quantum numbers of the subgroups and the Car- 
tan subalgebra. Hence, we have shown that the general 
Hamiltonian so that the general Laughlin states are the 
non-degenerate ground state. It is obvious to see that 
it preserves the conservation law of the center-of-mass 
position. 

So far there is no efficient way to calculate the Clcbsh- 
Gordan coefficients for the general SU(N)xSU(N) to 
SU(N) case. In the SU(3) case, only the Clebsh-Gordan 
coefficients for small multiplets can be calculated numer- 
ically. However, F(R, x) can be obtained analytically in 
general because 



(mp,0,0,..,0)iv x (mp,0,0,..,0)jv 

mp 

= 0(2mp-2g, g ,O,..,O) JV , 

9=1 



(42) 



where the right hand side contains no repeated multi- 
plets. We leave this important mathematical problem 
to mathematical physicists and conjecture that the pair- 
hopping range defined by the half width of the \F(R, x)\ 2 
is also long-ranged scaling as yfp with an exponential de- 
cay tail for the general case. 



IV. SINGLE MODE APPROXIMATION 

In the SU(2) case, Girvin et al. 10,11 showed that for 
any liquid ground state in the lowest Landau level, the 
density fluctuation costs finite energy, which implies in- 
comprcssibility, within the single mode approximation. 
In this section, we generalize their result to the SU(N) 
model. 

As mentioned in the earlier section, our lattice model 
corresponds to the fractional quantum Hall effect in 
QpN-i wrncri j s su bject to the background U(l) mag- 
netic field with the quantization n. The single-particle 
Lagrangian in the lowest Landau level is given by 



-inip a 



dlpg 

dt 



(43) 



Using Eq.(32), in the flat-space limit Eq.(43) can be writ- 
ten as 



N-l 

fc=i 



ny k x k = AjXj 



(44) 



where we set \z k \ -C 1 and z k = x k + iy k for fc = 1 to 
N — 1 and Aj and Xj are the 2(N— l)-dimensional vector 
potential and the position vector {{x kl y k }) respectively. 
From Eq. (44) , the single-particle orbit in the lowest Lan- 
dau level can be obtained as 



JV-l 



i 



k=l 



v / 2n2 l «:l k 



(45) 
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where l k are non-negative integers. We recognize that 
Eq.(45) is the product of the (N—l) lowest-Landau-level 
wavefunctions in two space dimensions. Thus, in the flat- 
space limit, the lowest Landau level in CP^ -1 becomes 
the direct product of (N — 1) quantum Hall effects in 
two dimensions, ft can be seen by the non-commutative 
algebra as well. From Eq.(44) 



/(fci, k 2 ) can be easily computed as 



[xk,Vk] = for k = 1 to N - 1 

n 



(46) 



which indicates the there are N—l independent non- 
commutative planes. In the rest of the section, we present 
the results for N = 3. The formalism can be generalized 
to any TV easily. We also note that in the rest of the 
section N is the symbol for the number of particles. 

In the single- mode approximation (SMA), the varia- 
tional excitation energy is given by 



A(fci,fc 2 ) = f(h,k 2 )/s{k 1 ,k 2 ), 



(47) 



where fcj are the complex wave numbers in the i th quan- 
tum Hall plane and f(ki,k 2 ) and s(k\,k 2 ) are the os- 
cillator strength and the static autocorrelation function 
given by 

f(ki,k 2 ) = — < * m \[p\ klMy [H,p (klM) ]]\^ m >(48) 
s{k 1 ,k 2 ) = ^ < ^m\p\ kuk2) P{ kl M)\^m > (49) 



respectively, where P( kl ,k 2 ) 1S * ne density operator. Here 
we use adopt Girvin et al.'s notation 10,11 which is a little 
bit different from our previous Letter 6 . In the lowest Lan- 
dau level, both f(k\,k 2 ) and s(k\,k 2 ) should be treated 
carefully because coordinates do not mutually commute. 
Particularly, the kinetic energy vanishes and they should 
be redefined by 

Kh,k 2 ) = i < ^ m \\p\ klMV [V,P {kuk2 )W m > (50) 
S{h,k 2 ) =-^<*m|p[ fclifca) P(fc 1 ,fc 2 )|*m> (51) 

where p is the projected density operator and V is the 
projected potential energy in the lowest Landau level, 
which are respectively given by 



P(ki,k 2 ) 



N 

£ 

72„ J2, 



g g 2 dz 2j g 2 Z 2} g 2 Zlj 



- 1 f d z qi d z q 2 f 
V= 2j ~(2n)^ vMiP l' 



bii 2 +b2i 2 , 
92 )P( gi ,«2)-pe 2 )(52) 



where p is the average density. In Eq.(52), 92) is re- 
quired to be positive indicating the repulsive interaction 
to ensure the excitation energy to be positive. Using the 
algebra for the density operator 

k*qi+k*q 2 k^+k^* 

[P{k u k 2 ),P(q u q 2 )] = ( e 5 - e 5 )P(k 1 ^ n ,k 2 +q 2 ), (53) 



_ 1 ^ — ^ q*k 1 +q*k 2 q 1 k*+q 2 k* 

f(ki,k 2 ) = - 2^ v{ qi ,q 2 ){e s -e 5 )[ 



2 ±\lc„\ 2 k*q-,+k 



fcll^ + |fc 2 



s(qi,q 2 )e 2 (e 



+s(fci+qi,fc 2 + g2 )(e 3 -e 5 )] (54) 

A direct expansion shows that j(k\,k 2 ) vanishes in the 
fourth order in k. To show the necessary condition for 
existence of the excitation gap, we have to demonstrate 
that s(k\,k 2 ) vanishes in the same order in k. Then, 
because of the isotropy between k\ and k 2 

v ; | 4 + d\kx\ 2 \k 2 \ 2 + c\k 2 \ i v ; 

remains finite as k approaches to zero in any direction, 
where a, b, c , and d are constants. 

The asymptotic behavior of s(ki,k 2 ) can be analyzed 
by relating with the radial distribution function g(r) by 



s(k) = 1 + pj d 4 re' 4fe - r [ 5 (r) - I ] + /9 (2 7 r) 4 5 4 (/c) (56) 

where k — (Refci, Imki, Rek 2 , lmk 2 ) and r — 
{x\,yi,x 2 ,y 2 ) are real vectors. For filling factor v = 
1/m 2 we obtained 



p\g{?) i] 



( 



< n h i 2 n 00 > - < n hh >< n o > -^ 2 -<^(; 1 ; 2 )(oo)) (57) 

where r\ = xj + y 2 and r\ = x\+y%, and n hh = cj^cj^ 
where cj j 1S the electron creation operator in the or- 
bit (h,l 2 ). Establish the relation between s(k\,k 2 ) and 
s{ki,k 2 ) 



s(k 1 ,k 2 ) = s(k 1 ,k 2 )-(l-e 2 ), 



(58) 



and define 

M nin2 ee J d 2 ri d 2 r 2 (|)"Hy)" 2 p[3(ri,r 2 ) - I] (59) 

One can compute easily that 

M 00 = m 2 (< Nn Q0 > - < N >< n 00 >) - 1 
Mio = m 2 (<(Li+N)n oa >-<L 1 +N><noo>)-l 
M i = m 2 (<(L 2 +i%oo>-<i2 + A r ><^oo>)-l(60) 

where N = j n hh ^ s the total number of particles and 
Li = X^ li2 Gi) n iii2 are the total angular momentum on 
the i th quantum Hall plane. Because of the conservation 
of the number of particles and the angular momentum, 
their fluctuation is zero, so M o = M w = M i = — I. 
Then, the second order of k in s(ki,k 2 ) vanishes. The 
asymptotic behavior of s(ki,k 2 ) indeed scales as the 
fourth order in k. Eq.(55) holds true. 
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The current analysis relies on the transformation from 
CP^ 1 to R 2 ^ -1 ), namely the flat-space limit. The for- 
mer one is a compact space whose volume is finite without 
boundary. The later one is a non-compact space whose 
volume is infinite but the exponential term in Eq.(45) 
sets the natural boundary. An importance question is 
whether the structure of the energy spectrum is pre- 
served in the transformation. We believe that the en- 
ergy spectrum is not a one-to-one mapping, because in 
p2(;v-i) t nere are certainly gapless edge excitations. On 
the other hand, in CP^ -1 there is no edge excitation due 
to the lack of the boundary. However, the structure of 
the bulk excitation is preserved because we do not intro- 
duce any flux which generates the electro-motive force to 
close/open a gap in the transformation. Therefore, the 
existence of the finite excitation gap in R 2 ^^ 1 ) implies 
that our lattice model also has a finite excitation gap. 
Besides the edge modes, whether or not there are other 
gapless bulk excitations can not be answered by the cur- 
rent approximation in R 2 ^ -1 ). If there are gapless bulk 
excitations in R 2 ^ -1 ), it should be also true in our lat- 
tice model. As far as we can say, within the single mode 
approximation, we conclude that our unique ground state 
describes an incompressible quantum liquid. 

V. DISCUSSION AND CONCLUSION 

Guided by the new symmetry of the center-of-mass po- 
sition, we construct a family of the models to describe the 
fractionally-filled incompressible liquid in any dimension. 
They are long range pair hopping model in the frustrated 
lattice in d > 2. We prove rigorously the uniqueness of 
the ground state in the open boundary condition using 
the group-theoretical method. We also compute the en- 
ergy gap using the single mode approximation, which 
is still the best analytical method to show the finite- 
ness of the energy gap in the fractional quantum Hall 
effect. Since the model is highly related to the higher di- 
mensional generalization of the quantum Hall effect, one 
can generalize our models in different topological struc- 
tures, for example, the torus and discuss the possibility 
of the fractionally-charged excitation for the future ex- 
plorations. 

We dedicate this work to Darwin Chang who made 
important contributions at the early stage of this work, 
but passed away before its completion. A special and 
deep gratitude should give to Dung-Hai Lee for leading 
the author to this topic. We are supported by the NSC 
97-2112-M-002-027-MY3 of Taiwan. 



Appendix A: Proof of the non- degeneracy of the 
zero-energy state 

In this section, we present the rigorous proof of the 
uniqueness of the ^ m as the zero-energy state of the 
Hamiltonian in the Eq.(37). The key procedure has been 



outlined in our previous Letter 6 . 

For simplicity, we focus on m = 3 and set p to be an 
odd integer. In this case, Eq.(37) has the following form 

H = — ^ p^p- 2 ' 1 ' ' '-' ) (Al) 

The Hamiltonian contains the term which projects two 
SU(N) spins in the {3p, 0, .., 0)n to the (6p — 2,1,0,.., 0)j\r 
subspace. The direct product of two SU(N) {3p, 0, .., 0)jv 
spins can be decomposed by the SU(N) subspace given in 
the Eq.(42). Particularly, when considering the antisym- 
metric combination, the complete set of spaces reduces 
to 

(mp, 0,0, ..,0) N x (rnp, 0,0, ..,0)jv| a 

mp 

= {2mp-2q,q,0,..,0) N , (A2) 

q— 1, odd 

where only odd q are allowed. If \ is t ne zero-energy 
state of Eq.(Al), there is no (6p — 2, 1, 0, .., 0) component 
in its two-spin spectrum. In other words, 

/p(6p-6,3,0,..,0) w + p(6p-10,5,0,..,0) N + + p(0,3p,0,..,0)x _ ^ 

(A3) 

for any pair {ij). Eq.(A3) constraints the symmetry 
properties in the zero-energy state. We note that the 
state in the SU(N) (3p, 0, .., 0)n multiplet can be writ- 
ten as the symmetric product of 3p SU(N) fundamental 
spinor in the Eq.(32). For particle j, it is given by 

( A4 ) 

where ctjk runs from 1 to N. Eq.(A4) is the alternative 
way of expressing Eq.(33). Using Eq.(A4) as the basis, \ 
in general can be written as 

JV d(p) 3p 

x= e cwnn^"' (as) 

= j=ln=l 

where the wavefunction C{{aj n }) satisfies the following 
Schrodinger equation from the Eq.(A3) 

C{..,{a l }..,{a J }, ..) 

= ]T L4 3 (K},K};{A},{&}) 

{ft}, {ft} 

+A 5 ({a i },{a j };{Pi},{Pj})+ - 

+A 3p ({ ai }, {«,}, {ft}, {&})]C(.., {Pi}, .., {&}, ..) 

(A6) 

for any pair {ij), where A q are the tensors for the pro- 
jection operator p(® p ~ 2q > q >°'--> 0> ) N . A q does the following 
symmetric operations 

1. q of the 3p indices of particle i is made antisym- 
metric to q indices of particle j. 
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2. the rest of the indices of particle i is made totally 
symmetric to the rest of the indices of particle j. 

There is degree of freedom to choose which pair of indices 
is made antisymmetric in the symmetry operation given 
above. For example, A 3 can be written as 

^(W.to); {&},{&}) 

= ^(^»C 1 -^C 1 )(C ,P+1 C' ,P+1 -C ,P+1 C' ,P+1 ) 

jVjj P' 1 ft' 1 Pi 1 P^ P«,P+i Pj,p+i Pj.p+i P»,P+i ' 

(C 2p+1 C ,2p+1 - C' 2p+1 C' 2p+1 )(C 2 -C p C p+2 -C ,2p 

v Pi, 2p+l Pi,2p + 1 Pj,2p + 1 Pi,2p + l'' v Pi2 Pip Pi, p + 2 Pi, 2p 



Pi,2p+2 Pi,3p 

+ sym.), 



J«i,3p^aj2 ja 3 p^ajp + 2 £aj,2p£aj,2p+2 ^ 

ft,3p ft'2" ftp ft, p+2" ft,2p ft,2p + 2 " i 



: a j,3p 
ft,3p 



(A7) 



where N 3 is the normalization constant. Consequently, 
for a particular pair (ij), one can arrange the antisym- 
metric pairs so that C becomes — C by the following in- 
dependent exchanges 



(aa..a ip ) «-» (aji..a jp ) 
(ai,2p+i--ai,3p) (aj,2p+i"Q!j,3p) 



(A8) 



If the symmetry property of C shown above can be made 
true simultaneously for any pair (ij). The argument that 
^ m is the unique zero-energy state can be given as the 
following. Let us consider the independent exchange of 
the first group of p indices while keeping others fixed. C 
becomes — C has to be established in any pair (ij). Be- 
cause the number of particle is N — d(p) is exactly equal 
to the total number of states that p indices represents, C 
is proportional to 



C ~ e{(a jl ..a jp )} 



(A9) 



where e{(aji..Oj P )} is the rank d(p) tensor with respect 
to the group exchange. Similarly, the second and third 
properties in the Eq.(A8) leads to 

C ~ e{(a jl ..a j p)}e{(a jtP+1 ..a jt 2p)}e{(a jt 2p+i..a jt3p )} 

(A10) 

Using Eq.(AlO), any zero-energy state \ is proportional 
to * m . 

Now, we shall prove that Eq.(A8) can indeed be made 
true for all pairs (ij) simultaneously. Let us assume that 
there exists a ground state solution whose C does not 
satisfy Eq.(A8) for pair (kl). It means there is at least 
one group exchange, say {oiki-.otkp} {ctn-.ai p }, so that 
C does not follow Eq.(A8). However, since the wavefunc- 
tion still has to satisfy Eq.(A6) for (k, I), one can write 



3p 



C= C q , 



(AH) 



q— 3, odd 



where C q is the component of C that is odd with respect 
to exchange of exactly q pair of indices between particle 
k and I and even with respect to the exchange of the rest. 
Now let us consider the effect of {aki-.oikp} {aii-ai P } 



on C. After the exchange, C q can either change sign 
or stay invariant depending on whether an odd or even 
number (out of q) antisymmetric indices are contained in 
the specified triplets. In other words upon {ctki-.oikp} <^ 
{an..aip} we have 



3p 
q— 3, odd 



(A12) 



where r\ q = ±1. Since Eq.(A8) is not satisfied, all r\ q 
must not simultaneously be —1. Now consider a new C 



c 



c 



3p 
q=3,odd 



(A13) 



It is obvious that upon {a k i..a kp } f-> {an..ai p } C — > 
—C. Moreover by construction C only contains those 
C q whose r] q = —1. Now use C as the starting C and 
repeat the above operation until we reach a final C for 
which Eq.(A8) holds for all triplet exchanges and for all 
(ij). Since at each stage of obtaining C certain C q are 
projected out, there must be missing components in the 
final C. However we have already proven that any C that 
satisfy Eq.(A8) for all (ij) pair must lead to the solution 
X ~ ^m- However, ty m contains all components for all 
pair (ij). Consequently we have reached a contradiction. 
Therefore it must be possible to make Eq.(A8) hold true 
for all pairs (ij) for any ground state solution satisfying 
Eq.(A3). 

The proof can be generalized to any m by assigning m 
groups of indices. Thus, we have proven that ^ m is the 
unique zero-energy state of Eq.(37). 



Appendix B: Summary of SU(3) algebra and 
representation theory 

1. Algebra 

The SU(3) group is the one with which people are very 
familiar besides SU(2). This note will not be a thorough 
review of Lie algebra but focuses on what we shall need 
in the paper. The generators of SU(3) 
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satisfy 



Tr{T a T h ) = -S ab 



(Bl) 



in the standard convention. The Cartan subalgebra con- 
tains T3 and T$. Denote them by Hi and if 2 respectively. 
The simply roots of SU(3) can be obtained as 



, 1 y/3 2 1 \/3 
« X = (2.-2-). « =(2'-^ } 



All positive roots of SU(3) are given by a 1 , a 2 , and a 1 + 
a 2 . The generators correspond to the positive roots are 
given by 



E, 



E a i 

E a 2 

a 1 +a 2 



1 
1 

71 
1 

= 71 



(T 4 + iT 5 ), 
(Te - iT 7 ), 



(B3) 



These generators are the raising operators in SU(3). 
Their lowing operators are the Hermitian conjugates of 
themselves, i.e. 



1 

71 



(T 4 - iT B ), 



E_ 



* 2 = y|( T 6 +«7V), 



V2 



The algebra is given by 

[E a i,E_ a i 



IE, 
[E, 



E 9 

a* 1 — d 



El — a 1 ■ H 
E 2 = a 2 -H 
1 



—^E a l +a 2 

[E- a i,E a i +a 2] = -^E a 2 
1 

[E-a 2 -, E a l +a 2\ — ——^E a 2 



The fundamental weight is defined by 
2a l ■ fi j . 



a 1 



In SU(3), /x"s are given by 



fj, 



V3, 



6 



Lb = -,- 

P V 6 



(B4) 

(B5) 
(B6) 

(B7) 
(B8) 
(B9) 

(BIO) 
(Bll) 



The representation whose highest weight is the funda- 
mental weight is called the fundamental representation. 
Since the rank of SU(N) group is N — 1. The number of 
simply root and that of the fundamental weight are also 
N — 1. The highest weight /U in any SU(3) representa- 
tion is given by // = put 1 + qji 2 . p and q are called Dynkin 



coefficients, which are unique for every representation. 
Therefore, SU(3) representations are denoted by (p,q). 
The total number of the Casimir operators of SU(N) is 
also equal to its rank. We shall pay our attention to the 
quadratic Casimir operator only. It is defined by 



C = Y J T a T a 



(B2) We can compute it to be 



C = H 2 + H 2 



all positive roots 



E n E- 



(B12) 



(B13) 



Using Eq.(B5)(B6)(B7)(B8)(B9) and the highest weight 
method, we compute the quadratic Casimir for SU(3) 
(p, q) representation: 

C\p, q] = \ (p 2 +pq + q 2 + 3p + 3q) (B14) 



The dimension of (p, q) representation is given by 

j% g ]= (p+i)(g+ y p+9+2) (bis) 



In the context, we introduce the SU(3) algebra in favor 
of particle physics. Namely, the SU(3) algebra is given 
in the Cell-Mann notation. There is another basis which 
is also very interesting and useful in certain problems 14 . 
If we define 



(B16) 



where i, k,/j,,is = 1, 2, 3, having the following properties: 



A\ 



(4 fe ) f 

: 



(B17) 
(B18) 



in which we know that there are only 8 independent gen- 
erators. It can be checked that they satisfy the following 
commutation relations: 



[AiAi] = siAi 



SlAl 



(B19) 



These 8 independent generators form the SU(3) algebra. 
The Cartan subalgebra is given by 



h 1 = \{A\-A 2 2 )=T, 
h 2 = 1 2 (A 2 -Al) = - 1 2 T 3 + ^T 8 



(B20) 



This notation is so-called Okubo's notation. The relation 
between Okubo notation and the Cell-Mann notation is 
given as the following. Denote 

T\ = h, T 2 = I2, T 3 = I 3 
T± = K Xl T 5 =K 2 
Ta = Li, T7 = L 2 

T$ = M (B21) 
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and 



I±=h± ih 
K ± =K 1 ± iK 2 
L± = L 1 ± iL 2 , 



(B22) 



and the A\ can be written as 



A\ = h + \>ftM, A\ = I+, A\ = /_, 

lV3M, Al=K+, A\ = K. 



A\ = -h 



A 3 3 = --VSM, Al = L + , Aj = L_ 



(B23) 



On the other hand, it is possible to obtain the transfor- 
mation between /, K , and L: Define 



Pk = e 



e 

i-kK 2 



(B24) 



Then 



p t ~ 1 i±p i = pr 1 K±p l = l±, 

P^I±P k = L T7 P^K ± P k = -K T , 
Pf l I±Pi = K^P^KiPt = -J±, 
P- X L ± P % = -K±, P^L ± P k = -I T , 

P^ l L±Pi = -L T (B25) 



2. Representation in the X-L 3 basis 

The X operator satisfies the following equation 

X\(p,q)jLL z ) = 2p + q-6j\(p,q)jLL z ) (B26) 

j ranges from 0, 5, ..to j — is the highest X state. 
j is a quantum number on the X-axis. However, due to 
degeneracy, for a certain j, it could be many L. The 
value of L is given by \j — §|...j + §. Using this basis the 
SU(3) Clebsch-Gordan coefficients (isofactor) are given 

by 

((p u Q)M,{p 2 ,0)N\\{p,q)jL) 

I (p + l)(q + l)!(p + q + L-M-N + l)l 



{-If 



(pi - ?) ! (P2 - q)l 



(p + q— L — M — N)\(2M + 2N - q + 1)! 
V { Pl -2M)\(p 2 -2N)\ 

V / (2M + 1)!(27V + 1)! 
(p + q+iy. 
V (p L -2I{)\(p 2 -2I' 2 )\ 



(2JJ)!(2J£)!(2M - 2/()!(2JV - 2I' 2 )\ 



1 1 1 2 2 

x { M-I[ N-I' 2 j 
M N L 



(B27) 



where j = M + N — | and we use the Wigner 9-j symbol. 
If L = M + N, namely L = j + | case, the Wigner 9-j 
symbol has a simpler form: 

1 i 1 2 2 

M-I[ N-I' 2 j 
M N M + N 



1 

~ V(l + !)( 2M + 1)(2^V + l)(2j + T] 
The CG coefficient becomes 

((p 1 ,0)M,(p 2 ,0)N\\(p,q)jL) 



(B28) 



= (-1) 



.. {p+l)g\{p + g-2M-2N)\ 



(pi - q)Kp2 - ?) ! 



/ (2M)\(2N)\(2M + 2N - q)\ 
(p + q + iy.( Pl -2M)!(p 2 -2N)\ 

V r ,x24 (PL- 2/()!(p 2 - 27 2 )! 



J l^ J 2 2 



(2I()!(2J£)!(2M - 2JJ)!(2JV - 2I' 2 )\ 

(B29) 



This result can be compared with the CG coefficient to 
the highest weight: 

((pi,0)M,(p 2 )N\\(p,q)0,^) 



(P + !)!«! 



1 ' (p + g + l)!(pi-g)!(p2-«)! 
(pi - 2M)!(p2 - 2JV)! 



(2M)!(2iV)! 



?M+N, § 



(B30) 
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